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SUMMARY 


A study is made of the advantages that can be realized in 
compressible-flow research by employing a substitute heavy gas in 
place of air. Heavy gases considered in previous investigations are 
either toxic, chemically active, or (as in the case of the Freons) have 
a ratio of specific heats greatly different from air. The present report 
is based on the idea that by properly mixing a heavy monatomic gas with 
a suitable heavy polyatomic gas, it is possible to obtain a heavy gas 
mixture which has the correct ratio of specific heats and which is non- 
toxic, nonflammable, thermally stable, chemically inert, and comprised 
of commercially available components . 

Calculations were made of wind-tunnel' characteristics for 63 gas 
pairs comprising 21 different polyatomic gases properly mixed with each 
of three monatomic gases (argon, krypton, and xenon). For a given Mach 
number, Reynolds number, and tunnel pressure, a gas -mixture wind tunnel 
having the same specific-heat ratio as air would be appreciably smaller 
and would require substantially less power than the corresponding air 
wind tunnel. Very roughly the results are as follows: 


Mixtures 


Power required relative 

to air tunnel . 

Argon 

Krypton 

Xenon 

0.3 

0.2 

0.1 

Characteristic dimension of 

test section relative to air tunnel . . 

0.7 

0.6 

0.5 


Analogous advantages can be realized in compressor research and in 
firing-range research. Other applications arise through the possibility 
of selecting and proportioning a gas mixture so as to have at wind-tunnel 
temperatures approximately the same dimensionless characteristics as air 
at flight temperatures. Characteristics which involve relaxation time 
and variation of specific heat with temperature fall within this category. 

Most of the polyatomic gases considered are more inert and stable 
than Freon 12. One of the more promising gases is bromotrifluoromethane 
(CBrFg), commercially called Freon, I 3 BI or Kulene I 3 I. Of the three 
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monatomic gases considered, all are completely inert, although only argon 
now is available in large quantities. 


INTRODUCTION 


The reasons for considering gases other than air as possible test 
media for compressible-flow research stem primarily from the relatively 
low speed of sound in certain gases. In general, the heavier the gas, 
the lower the speed of sound at a fixed temperature. Hence, experiments 
conducted at a given Mach number in a heavy gas will be conducted at lower 
velocity than the corresponding experiment in air. Some significant 
advantages of conducting wind-tunnel, firing— range, and compressor exper- 
iments at a reduced velocity have been noted previously in the work of 
Theodorsen and Regier (ref. l) , Smelt (ref. 2), Kantrowitz (ref. 3), 

Huber (ref. 4 ), Buell (ref. 5), Donaldson (ref. 6), and von Doenhoff and 
Braslow (ref. 7). Some uses of gases other than air in shock tube research 
are discussed by Duff in reference 8. These investigations show that the 
substitution of a heavy gas for air offers the possibility of: (l) extend- 

ing the range of existing research apparatus; (2) achieving greater econ- 
omy of construction and operation of large high-speed wind tunnels; and 
( 3 ) providing greater facility in obtaining data for special types of 
research. These advantages will be discussed briefly in the order listed. 

The possibility of extending the range of existing apparatus was 
clearly demonstrated by the experiments of Theodorsen and Regier, which 
appear to be the first compressible-flow experiments in which supersonic 
aerodynamic data were obtaijied in a gas other than air. By rotating pro- 
pellers in Freon 113 (CCI2FCCIF2, having a speed of sound 0.39 times that 
of air, and a density 6.5 times as great) they were able to achieve tip 
Mach niunbers of 2.7; whereas the highest tip Mach number achieved with 
the same apparatus using air was 1.0. Also, much higher Reynolds numbers 
were obtained. In a similar fashion, Buell employed the Freons, and 
Donaldson employed xenon as the test medium through which projectiles 
were fired in order to obtain Mach numbers much higher than could be 
obtained with the same apparatus by firing through air. An. extended range 
of operation of a wind tunnel has been demonstrated by von Doenhoff and 
Braslow who report that the maximum attainable test-section Mach number 
of the Langley low-turbulence pressure tunnel was increased from 0.4 to 
1.2 by replacing air with Freon 12. 

A possibility of achieving greater economy by employing a substitute 
heavy gas in a wind tunnel arises because the test-section dimensions are 
smaller and the horsepower required is only a fraction of that required 
for an air wind tunnel operating at the same Mach number, Reynolds number, 
and pressure. Smelt has made a general study pointing out this possibility 
for various inert gases that have, unfortunately, a specific-heat ratio 
different from air. The importance of reduced power requirements, if 
accomplished without sacrifice in ratio of specific heats, needs no 
elaboration in view of the power of modern wind tunnels. 
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The third advantage mentioned, namely, facilitating certain types of 
aerodynamic research, stems directly from the lower test velocity relative 
to air. One consequence of reducing the velocity while maintaining physi- 
cal similarity of a flow is that centrifugal stresses of rotating machin- 
ery are reduced as the square of the velocity. This can be an important 
advantage in the design of large wind-tunnel compressors. Also, for a 
given Reynolds number and Mach number, the stream dynamic pressure is 
significantly reduced. Consequently, it is much easier to construct cer- 
tain models for testing in a heavy gas than for testing in air. Consid- 
erable compressor research has already been conducted in Freon 12 which 
would not have been practical if air were employed as the test medium. 

(See the work of Kantrowitz for example) . The feature of a reduced 
dynamic pressure also would facilitate tests in pressurized transonic 
wind tunnels wherein model stresses or model-support stresses limit 
undesirably the useful angle-of -attack range. Firing -range research is 
facilitated through use of a heavy gas since a lower launching velocity 
results in longer time intervals for observation, as well as lower model 
stresses during launching. Flutter research may be facilitated because 
lower velocities result in lower oscillation frequencies, which, in turn, 
require simpler instrumentation for certain measurements. Wind-tunnel 
bomb-drop tests are facilitated because the lower velocity enables a 
smaller model to be used while simulating the required value of the Froude 
number (uoo^/gZ) . For example, simulation of the Froude number at a Mach 
number of one requires a model practically as large as the prototype if 
air is employed, but only about l/4 scale if a gas four times as dense 
is employed. 

As regards the various disadvantages involved by using substitute 
heavy gases in aerodynamic research, it is clear that even without detailed 
study these are more obvious than the advantages. The main disadvantages 
are: (l) The ratio of specific heat (7) differs from that of air for 

the heavy gases proposed thus far] (2) the use of any gas other than air 
results in more severe practical problems of operating a research ..facility; 
and (3) conversion of existing apparatus for operation with a heavy gas 
can be costly.^ It is to be noted that gases with 7 significantly 
different from air would not yield in a supersonic flow at a given Mach 
and Reynolds number, either the same skin friction or the same pressure 
distribution as air would yield. Further elaboration here would not be 
pertinent inasmuch as the present report concerns a method of eliminating 
the first and chief disadvantage while retaining the advantages previously 
discussed. 

^An additional disadvantage sometimes expressed (e.g., ref. .9) is 
that the relatively long relaxation time for molecular vibrations in 
certain polyatomic gases would invalidate data obtained at supersonic 
velocities. However, in the case of most gases considered here, this is 
believed not to be the case. Reasons for this are based on calculations 
presented in an appendix and are discussed later. 
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There are many rather evident practical disadvantages of operating 
a facility with a gas other than air. Purging of the test-section region 
(or entire wind tunnel) before each model change handicaps the operation 
of any research apparatus whether it is a wind tunnel, compressor test 
facility, or firing range. T?he necessity of preventing excessive contam- 
ination by air and of minimizing leakage of the heavy gas to the atmos- 
phere is a definite practical disadvantage. Also the operation of the 
necessary purification equipment, the increased pumping equipment, arid 
the additional gas storage equipment represent such disadvantages. 

As for the third disadvantage, namely, that involving the cost of 
converting existing equipment for operation with a gas other than air, it 
would appear an insurmountable one for tunnels that are not pressurized, 
and it would be a dominant consideration for existing large wind tunnels 
in which the compressor speed cannot be varied efficiently. To maintain 
proper matching of compressor and tunnel, the lower stream velocity of a 
heavy gas would require that the compressor be operated at correspondingly 
lower speed. Alterations to accomplish this can be costly, as can be the 
possible additional equipment required for pumping, purifying, and storing 
the heavy gas. 

A general appraisement of the various advantages relative to dis- 
advantages cannot be made since this would depend on each particular 
case. The purposes of this report are to indicate a method of construct- 
ing a satisfactory heavy gas for aerodynamic testing and to evaluate the 
magnitude of the advantages such gases offer. Knowledge of these advan- 
tages will enable a comparison with the disadvantages to be made in a 
given case. 

The present research was begun upon conception of the following 
simple idea: The proper ratio of specific heats (7 = 1 . 4 ) can be achieved 

with a heavy gas by suitably mixing a heavy monatomic gas (7 >1.4) with 
a heavy polyatomic gas ( 7 < 1 . 4 ). Exploitation of this simple idea turned 
out to be unexpectedly arduous because of difficulties involved in obtain- 
ing adequate chemical, physical, pharmacological, and thermodynamic data 
on the many polyatomic gases known to modern chemistry. Primarily as a 
result of recent developments in the field of fluorochemicals, it was 
possible to find among the known polyatomic gases over 30 which were 
indicated by available data to be sufficiently nontoxic, nonflammable, 
chemically inert, and thermally stable to justify consideration. Evalua- 
tion of the relative advantages of each polyatomic gas for aerodynamic 
testing requires calculating the thermodynamic properties from spectro- 
scopic data, the proper proportions for mixing from thermodynamic proper- 
ties, and then computing the density and viscosity that would result when 
mixed with each of three monatomic gases considered (argon, krypton, and 
xenon) . The remainder of this report describes the detailed methods 
employed in the evaluation, the results obtained, and the conclusions 
drawn. 
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After completion of the main body of calculations, a short unpublished 
note by E. F. Relf, in England, became available in which the possibility 
of obtaining the proper ratio of specific heats through mixing gases also 
was pointed out. Relf made a rough estimation for a mixture of sulfur 
hexafluoride and xenon, which is one of the 63 mixtures considered in the 
present research. 


METHODS OF SELECTING GASES AND COMPUTING 
CHARACTERISTICS OF VARIOUS GAS MIXTURES 

Selection of Gases 


The selection of monatomic gases for consideration is simple, since 
only four such gases are known - argon, krypton, xenon, and radon - that 
have a molecular weight greater than air.. Radon, the heaviest of all, 
unfortunately must be excluded because of its radioactivity. The remain- 
ing three gases are completely inert. Some of their properties pertinent 
to the present analysis are listed in table I. Although argon is the only 
one now commercially available in large quantities, both krypton and 
xenon will be considered since they might be of use for special research 
in small apparatus, and since their availability undoubtedly will increase 
in the futixre.^ 

The selection of polyatomic gases for consideration was rather 
tedious. Of the approximately 13,000 polyatomic compounds listed in the 
1952 Handbook of Chemistry and Physics (ref. 12), 103 were found to be 
gases at 0° C. Compounds that are not gaseous at 0° C were not considered. 
Most of the gases were not considered further because they possessed some 
intolerable characteristic such as being highly poisonous (e.g., COCI2, 
^2^10^ C2N2, SO2) , toxic (e.g., CH3F, CF2Br2, CH3CH2F) , combustible 
(e.g., COS, all hydrocarbon gases), unstable (e.g., CHCCl, CIN3, PSF3), 
or generally active chemically (e.g., BF3, F2O, CF3SF3, CF3OF) . Still 
others were eliminated for less obvious reasons, such as reacting with 
water vapor to yield corrosive products (e.g.,SiF4, PF3, COFg, TeFe), or 
having a tendency to polymerize under pressure (e.g., CCIFCF2, CH2CF2, 
CF2CF2)^ or being active when in the presence of mercury (F3CNNCF3), or 
imparting, upon inhalation, an offensive odor to the breath which persists 
long after inhalation (TeFg) • This elimination process left I6 polyatomic 
gases for further consideration. It was noticed that all but two (CO2 
and N2O, the lightest two on the list) contained at least two atoms of 
fluorine per molecule. Because fluorine appeared to be the key to the 

^Argon now costs about 10 cents per cubic foot of free gas in small 
quantities, and about 7 cents in large quantities. Because very large 
amounts of argon are required in the production and fabrication of tita- 
nium, the cost of argon in the near future may be greatly reduced (see 
e.g., ref. 10). Krypton and xenon now are produced in the United States 
only as by-products of air liquifaction, but recently have been produced 
on a more economical scale in Europe from special krypton plants ( see 
ref. 11) . 
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problem, recent technical publications on fluorine chemistry were scanned 
(especially refs. 13, 1^, and I 5 ) for further information. From this 
search 15 additional polyatomic gases were found which met the required 
conditions of being nontoxic, noncombustible, chemically inactive, and 
thermally stable. Among these 15 recently synthesized gases were several 
(e.g., CBrF3) which turned out to be the most efficient for wind-tunnel 
use. Undoubtedly, additional new polyatomic gases will continue to be 
synthesized in the future inasmuch as a large industrial effort is being 
expended on fluorine chemistry. Quite a few apparently satisfactory gases 
(e.g., C3CIF7 and C2HF5) are known to be omitted from the list. Conse- 
quently, the present report which primarily indicates current possibil- 
ities of gas mixtures in reducing wind-tunnel power is by no means 
exhaustive and might considerably underestimate possibilities existing 
several years hence. 

It is pertinent to focus attention here on the extreme properties of 
fluorochemicals . Their chemistry often is totally unlike that of analo- 
gous halide compounds. Fluorochemicals include some of the most toxic 
compounds known, yet also include the most inactive and nontoxic compounds 
known. For a given molecular weight, they have the lowest boiling points 
of known compounds, being paralleled in this respect only by the inert 
monatomic gases. Intermolecular forces and surface tension are extremely 
low. The thermal stability and resistance to electric breakdown of cer- 
tain fluorochemicals are, extremely high. Most of these extreme properties 
can be illustrated by a few examples. First consider the following two 
gases: chloroform, CHCI3, and fluoroform, CHF3. The toxicity of chloro- 

form is well known. Fluoroform, however, is indicated to be completely 
inert physiologically (ref. I6) . Guinea pigs can live in an atmosphere 
of 80 percent fluoroform and 20 percent oxygen without being affected 
either during the test or afterward. Sulfur hexafluoride, SFe, also is 
indicated to be physiologically inert on the basis of similar tests with 
small animals (ref. 17)* On the other hand, the deceptively similar 
sulfur pentafluoride, (SF5)2^ is extremely poisonous, since, for concen- 
trations as low as one part per million - in which concentration phosgene 
is harmless - sulfur pentafluoride is lethal (ref. I8) . Further illustra- 
tion of the extreme properties of fluorochemicals is provided by carbon 
tetrafluoride, CF4. The thermal decomposition temperature of CF4 has not 
yet accurately been determined since it is so high as to be in the vicinity 
of the carbon arc temperature (ref. I9) • Carbon tetrafluoride has been 
heated in the presence of many metals, including Cu, Ag, Hg, C, Si, Ni, 

Sn, Pb, P, As, Sb, Cr, Mo, W, S, Se, Fe, Ru, and Pt up to the temperature 
at which glass softens without any reaction occurring. (See ref. lU, 

P* ^ 33 •) Sulfur hexafluoride is equally remarkable in its chemical inac- 
tivity at all temperatures below the softening point of glass. 

The chemical symbol, molecular weight, chemical name, and trade name 
of 31 polyatomic gases considered as possible components for aerodynamic 
test media are listed in the first four columns of table II. A symbol 
designating the general classification of each gas is listed in the fifth 
column. The meaning of each symbol is explained in the footnote of the 
table. Listed in table III are certain data, to be utilized later. 
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regarding the boiling point, critical constants, viscosity coefficient, 
and spectroscopic frequencies of the various polyatomic gases. The 
soiorces of data for each gas are indicated by reference numbers in the 
last three columns of table II. These data were compiled from references 
13 to 16 and 18 to 56. 


Computation of Thermodynamic Properties 


first step in evaluating aerodynamic qualities of any gas is to 
determine its thermodynamic properties such as specific heat and enthalpy. 
When two gases are mixed, the resulting thermodynamic properties can be 
calculated from the thermodynamic properties of the individual components. 

subscript 1 to denote the monatomic component, and subscript 2 to 
denote the polyatomic component, the following equations apply for binary 
mixtures of thermally perfect gases: 


molecular weight 

m = ximi + X2m2 

(1) 

equation of state 

p = PRT/m 

(2) 

specific heat 

Cv + R = Cp = xiCpj_ + XgCpg 

(3) 

ratio of specific 

Vi£^q4-<-» 

Xi(Cp^ - R) + X2(Cp^ - R) 

(4) 

enthalpy 

H = XiHi + X2H2 

(5) 


where R is the universal gas constant (I.987 cal. mole"^), and x 

denotes the mole fraction of a component (equal to the fraction by volume, 
or the fractional partial pressure of the component) and is subject to the 
relation xi + x^ = 1. For all monatomic gases. 


Cpi “ " 2 ^ ~ cal.°K"^ mole"^ 


Hence, from equation (It-), it follows that for 7 = 7/5, 



( 6 ) 
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The thermodynamic properties of a polyatomic gas can be accurately computed 
from infrared and/or Raman spectroscopic data, using equations developed 
by methods of statistical mechanics (see ref. 57 ; for example). For 
linear molecules (e.g., CO2 and N2O) , wherein the atoms are arranged in a 
straight line. 



( 7 ) 


3 n -5 



1 


hcwi 


kT 



1 


( 8 ) 


The various wave numbers Wi determined from spectroscopic data for each 
gas are listed in table III. The constant hc/k is equal to I.4385 cm 
For nonlinear molecules not subject to internal rotation (e.g., CF^, CHFg, 
SFg , etc . ) 


^P2 

R 


4 + 




sinh^ 



3 n -6 



1 


hccji 


kT 



( 9 ) 


(10) 


For nonlinear molecules subject to a torsional mode of oscillation (e.g., 
CF3CF3, CHF2CCIF2; and all other ethane-like molecules) the above sum- 
mations were employed only for the 3n-5 normal modes, and the torsional 
contribution was computed from the Pitzer-Gwinn tables (reproduced in 
ref. 58) using values for the potential barrier estimated from spectro- 
scopic data. Computations were greatly facilitated through use of the 
tables of reference 58. Typical curves showing the molar specific heat 
as a function of temperature for several of the various gases considered 
are presented in figure 1 . Curves of this type, together with equation 
(4), were used for each gas to determine the mixture proportions necessary 
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to ottain a value of 1.4 for 7. Similar data for the enthalpy were used 
to compute Mach number, pressure ratio, etc., for the different flow con- 
ditions considered. 

It is noted that the thermodynamic properties calculated from spec- 
troscopic data agree well with direct calorimetic measurements. This is 
illustrated in the following table: 

Calculated Cp Experimental Cp 

(from spectroscopic (from direct calori- 
Gas T,°C data of table III) metric measurements) 


C02 

-75 

7.71 

7.78 ^ 

1 



20 

8.82 

8.80 

r ref. 

57 


220 

10.62 

10.63 

1 


CClgFg 

-30 

15.6 

15-7 " 



0 

16.6 

16.7 


51 


45 

17.8 

18.0 1 

^ ref. 


90 

18.9 

19-1 J 



CF3CF3 

-50 

21.1 

21.0 

1 



10 

24.7 

24.4 

y ref. 

22 


90 

28.7 

28.2 J 

1 



Computations of Viscosity 


It will be seen subsequently that for a given Reynolds number and 
Mach number, the power required by a wind tunnel is directly dependent 
on the gas viscosity coefficient. The viscosity of a mixture of two gases 
can be approximately calculated from the equation (see ref. 59): 


1 ^ = 



+ 


1 + 



^2 



4 / maV/^ 


( 11 ) 


It is noted that a simple linear interpolation between the values of n 
for each component is not an adequate method of estimating the viscosity 
of the mixture. This is illustrated by several examples in figure 2 . 
Curves of this type were constructed for all gas mixtures considered. 

Experimental data are available for the viscosity of the three 

monatomic gases, but very little data are available for the various poly- 
atomic gases. Fortunately, the coefficient of viscosity can be estimated 
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quite satisfactorily from a knowledge of the critical temperature 
critical volume Vg, and molecular weight by the equation 


= 




(Vc)' 


2/3 



( 12 ) 


This equation may be deduced from the general Chapman-Cowling equations 
for gases having intermolecular forces dependent only on powers of the 
distance between molecules. (See refs. 59 3-nd 60.) Rather than to use 
the elaborate theoretical evaluation of the function F(T/Tc)^ the more 
simple empirical evaluation of reference 61 is used: 


F 



r 27.1 T/Tc 

27.1 T/Tq 

1 + 0.0952 Q- - 1.6^ 


for T<'1.6 Tc 

(12a) 

for T > 1 . 6 Tc 


where* T is in m is in gm./mole, Vq is in cc/gm. mole, and |j. is 
in micropoise. A comparison of the calculated and experimental viscosity 
coefficients of all gases considered for which experimental data could be 
found is given in the following table: 

Viscosity at 0^ C in micropoises 


Gas 

Calculated 

Experimental 

Ar 

210 

210 

Kr 

229 

230 

Xe 

211 

210 

SFq 

IU6 

142 

CCI2F2 

Il4 

118 

CHCIF2 

119 

120 

CO2 

134 

137 

N2O 

132 

135 

Air 

170 

172 


This agreement is satisfactory. Consequently, in the calculations of 
wind-tunnel power requirements, experimental values of |i are employed 
for the 9 gases listed above, but values computed from equation (12) are 
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employed for the remaining gases for which experimental data are not 
available. Since molecules of the above 9 gases are either simple or 
approximately spherical in structure, the accuracy of calculations from 
equation (12) is expected to be less for greatly different molecular 
structures, such as that of C4 Fio> but probably is sufficient for present 
purposes . 


Wind-Tunnel Power 


The general equation for compressor power required to maintain a 
steady rate of mass flow w of a perfect gas through a compression ratio 
CR starting with an initial temperature is giveri by the equation 


HP = J 


wc- 




(cr) 


{y-x)/y 


- 1 


(13) 


where J is a numerical constant depending only on the units employed, 

Tj is the combined adiabatic and mechanical efficiency (assianed to be 
independent of the gas), and Cp is the specific heat per unit mass. It 
is noted that Cp in equation (13) is related to the molar specific heat 
Cp and the molar universal gas constant R through the equation 



(Ik) 


By using a bar to denote a quantity divided by the corresponding quantity 
for air (e.g., m=m/(m)aij., Iff s HP/(HP) etc.), there 

results 


■HP = I TtY(7,CR) = T^Y(7,CR) 


(15) 


where Y(r,CR) s [^(cr) 


(r-i)/7 


, and L is a characteristic 


dimension of the test section. A more convenient equation is obtained by 
introducing the Reynolds number 
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HP = TtY(7,CR) (16) 

m 


from which it is apparent that,, as Smelt (ref. 2) has pointed out, power 
economy can he achieved by employing a gas with low viscosity and high 
molecular weight, or by testing at the lowest temperature possible, or 
by testing in the smallest facility (highest operating pressure) that 
will yield the given Reynolds number. The quantities Ke and L are 
related to the relative operating tunnel pressure p^ 



or, finally. 


_ Re|i 

L = Z 

Pt -/S 


(17) 


where 


f(7,Mj 


^^PcoUco 

PtH 


■) 


Equation (I7) shows that by increasing the molecular weight, or decreasing 
the viscosity, the size of wind tunnel required to obtain a given Reynolds 
number can be reduced. Hence, an alternate equation for relative wind- 
tunnel power is 


m ^ Y(7,CR) 


. (18) 


Equations (I6), (1?) and (18) are a generalization of similar equations 
developed by Smelt. 



NACA TN 3226 


13 


For subsequent computations, the above equations can be simplified. 

It will be assimed that CR depends only on (To achieve this in 

ah air tunnel, the compressor speed would have to be reduced for a heavy 
gas, thereby introducing the practical problem mentioned earlier.) For 
the case of transonic wind tunnels CR is the__order of I.3 or less, and 
direct numerical calculations_show that 0.97< Y(7,1.3)< 1 for any 7 
between 1.1 and l.U. Hence, Y = 1 is a good approximation for transonic 
wind tunnels irrespective of 7* For the case of supersonic wind tunnels, 
only gas mixtures having 7^i:l.lj- are considered, for which Y(1.4,CR) = 1. 
Consequently, in all cases of concern here, Y = 1 _is an adequate approxi- 
mation. A similar argument shows that the ratio. Y(7,CR)/ f(7,Mj is 
approximately unity for the range of 7 and considered here . Inasmuch 
as all calculations will be made on the basis that = 1, the above 
equations simplify to 


m 

(l6a) 

L 

Pt^/S 

(17a) 

: (^)^(IT)2 

(I8a) 

p^(m)®/2 


Four illustrative cases will be considered; (l) given Men, Re, and L; 

(2) given L, and p^; (3) given L, and HP; and (4) given mJ 

Re, and p^.__ Since for all cases Tt = M«, = 1, the reduction in velocity 
is always u«, = (m)”'*''^ . 

Case (1); Given Mo Re, and L .- This case corresponds to operating 

a given wind tunnel (L = l) with a gas other than air at the same values 
of and Re as for air. The horsepower required for the gas mixture 
relative to that for air is HP = ]l/m. However, from equation .( 17a) , 
it is seen that p.^ = for given values of Re and L. Hence, in 

order to achieve the same Re with a heavy gas as with air, the given 
tunnel would be operated at a lower pressure, resulting in model loads 
that are proportionately lower, since q = p.^. = IT/z/s. 

Case (2): Given L, and p-f" This case corresponds to substi- 

tuting a gas other than air in a given wind tunnel and operating at the 
same pressure ^d Mach number as for air. From equation (17a) it follows, 
that Re = ^/m/ This indicates that higher Reynolds numbers would be 
achieved with a heavy gas than with air, while (according to eq. (l6a)) 
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simultaneously drawing less power by the ratio HP = Since 

= 1 , model stresses are unchanged for this case. 

Case ( 3 ); Given tU, L, and HP .- This case corresponds to replacing 
air in a given tunnel by a gas mixture and drawing the same power for a 
fixed Mach number. The Reynolds jiumber , from equation (l 6 a), would be 
increased by the factor Re = m/h. In order to draw the same power as 
with air, the pressure and, hence, model loads would have to be increased 
by the factor p^ =//S. 

Case ( 4 ); Given M«a Re, and pt .- This case corresponds to conditions 
that might be prescribed in the initial design of a new wind tunnel whose 
size is to be determined by the required values of Mo, and Re, by the 
selected value of p^, and by the gas employed. From equation (l 8 a) it 
is evident that in this case, the horsepower required for a heavy gas 
mixture is less than for air by the factor 


HP 



(19) 


This represents a greater power reduction than in the cases above because 
the gas-mixture wind tunnel, according to equation (17b) , would be smaller 
than the corresponding air tunnel by the factor 


L = n/'/im 


( 20 ) 


The gas -mixture wind tunnel would produce the same aerodynamic data as 
the larger air tunnel with the same model stresses. Computations based 
on these latter two equations are presented later. 


Mach Number Limit for Condensation 


Although any of the polyatomic gases would be usable in apparatus 
where low temperatures are not involved, such as in a subsonic wind tunnel, 
a firing range, or a compressor research apparatus, ^ only a limited number 
would be useful in a supersonic wind tunnel where the static temperatures 
encountered are low and can result in condensation. An essential step in 
evaluating various polyatomic gases, therefore, is to determine the 
approximate useful Mach number range for each gas. The method of deters 
mining this was to draw the isentropic expansion curve for the partial 
pressure of each polyatomic gas on log-log paper together with the curve 
for saturation vapor pressure of that gas. As indicated in the following 
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sketch, the intersection point yields T^i n^ from which the maximum Mach 



number can be calculated from the relation 



( 21 ) ■ 


where 7 and h/T are determined from spectroscopic data by equations (4) 
through (10). Graphical solutions of this type were made for operating 
conditions of approximately 40° C (100° F) total temperature and several 
atmospheres total pressure. These solutions enabled the appropriate 
polyatomic gases to be selected for each of several design Mach numbers 
considered in the evaluation. It is emphasized that the graphical method 
included consideration of the variation of specific heat with temperature. 
As will be seen subsequently, such variation is of dominant importance 
in determining the limiting Mach number to which certain gases can be 
expanded without reaching saturation. Simplified criteria, such as boil- 
ing temperature,' are quite inadequate for determining the saturation limit. 
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RESULTS AM) DISCUSSION 
Wind-Tunnel Power and Size 


Complete computations could not be made for 10 of the 31 polyatomic 
gases listed in table II, since the necessary spectroscopic data and 
critical constants were not available for these 10 gases. Calculations 
have been made, however, of the relative wind-tunnel power requirements 
for 63 gas pairs comprising mixtures, of the 21 polyatomic gases for which 
data are available with each of the 3 monatomic gases. Mixture propor- 
tions were determined from equation (4) using the requirement that 
7 ^ = i.lf for flow in the test section. The calculations were made for 
three design Mach nxmbers, I. 3 , 2.5, and 3*5* Results are tabulated in 
table IV(a), IV(b), and IV(c), respectively. It should be noted that 
all 21 polyatomic gases are included in table IV(a), since all are free 
from condensation at M^, = 1.3* Only 7 gases for Moo = 2.5 (table IV(b)) 
and 2 gases for M„ = 3-5 (table IV(c)) have sufficiently low boiling 
points to avoid saturation for the assumed condition of 40 C total tem- 
perature. Operation at Mach numbers above about 3*5 would require total 
temperatures higher than the value of UO° C arbitrarily assumed. 

The values of relative horsepower HP given in table IV represent 
case (4) mentioned previously. This case corresponds to a comparison 
with an air wind tunnel for the same M^ Re, and p-^. If desired, any 
of the other three cases mentioned, or any of the individual relative 
quantities such as u„ or can be readily calculated from the values 
of m and ^ listed in the tables and from the appropriate equations 
developed earlier. 

The tabulated results show that there is no single polyatomic gas 
which is best for all Mach numbers or for use with all monatomic gases. 

For each monatomic gas, though, there are several polyatomic gases which 
appear roughly equal in their ability to require low wind-tunnel power. 

Of the gases that now are commercially available and not classified ques- 
tionable in table II, the most efficient ones for use with argon are: 

CH 2 F 2 , CBrFg, and CCI 2 F 2 for M« = 1.3; CBrFg and SFg for M„ = 2.5; and 
CF 4 for M«, = 3 . 5 . In general, the horsepower and size of gas-mixture 
wind tunnels relative to equivalent air wind tunnels are very approximately 
as follows: 


Mixtures 



Argon 

Krypton 

Xenon 

Power relative to air> HP . . . 

0.3 

0.2 

0.1 

Test-section dimensj^on 

relative to air^ L 

0.7 

0.6 

0.5 


These figures imply significant advantages of gas mixtures as wind-tunnel 
test media. 
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It is interesting to note that the heaviest polyatomic gas - C4F10, 
having = 8.2 - does not result in either the heaviest mixture or the 
best mixture. This is due to the high molar specific heat of which 

requires that only a small amount (6.5 percent) be mixed with a monatomic 
gas in order to obtain 7 = 1 . 4 . However, for applications to low- speed 
subsonic research where the value of 7 does not matter and, hence, a 
mixture is not needed, the heaviest gas is the best. For example, pure 
C4 Fio^ having 7 = I.06, would require only 1.8 percent of the power, and 
an apparatus one-fourth the size of that required by air for a given 
Re, and p-j^. In this respect pure C4F10 is several times more efficient 
than pure Freon 12 . (See table III where corresponding values are given 
for all pure polyatomic gases considered.) It is evident that pure poly- 
atomic gases having 7^tl.l are more efficient in reducing wind-tunnel 
power requirements than are gas mixtures having 7 = 1 . 4 . 

A possibility that should not be overlooked is that a value of 7 
close to 1.4 may not be necessary for transonic wind tunnels or other 
transonic research apparatus. In view of the surprisingly small differ- 
ences observed at ^transonic speeds by von Doenhoff and Braslow (ref. 7 ) 
between .uncorrected results in Freon 12 (7 = 1 . 13 ) and in air (7 = 1 . 4 ), 
it would appear that a reduced value of 7, perhaps between I.3 and 1 ^ 2 , 
might yield transonic data directly applicable to air for practical pur- 
poses. If a considerably reduced value of 7 is satisfactory in tran- 
sonic wind tunnels (such is not anticipated in supersonic wind tunnels), 
then the possible advantages of gas mixtures for transonic speeds are 
greater than the above calculations for 7 = 1.4 would indicate. The 
following table for CBrFs-Ar mixtures at Moo = 1.3 illustrates this: 



Case ( 4 ) 
Re, pt. 

: given 

and Moo 

Case ( 3 ) 5 given 
HP, Moo, and L 


7 

iff 

L 



TJi 

0.37 

577)+ 

2.0 

(from table IV(a)) 

1.35 

.30 

.69 

2.3 


1.3 

.22 

.62 

2.8 


1.25 

.16 

.53 

3.4 


1.2 

.10 

.44 



1.15 

.06 

.37 

6.2 

(pure CBrFg) 


These figures, which are typical of many of the argon mixtures, show that 
by reducing 7 to about 1.28, the power requirement for an argon gas 
mixture at transonic speeds would be about one-half that for the same type 
mixture proportioned to yield a 7 of 1 . 4 . This corresponds to one-fifth 
the power of an equivalent air wind tunnel operating at the same Reynolds 
number, total pressure, and Mach number (case ( 4 )).. Also, the Reynolds 
number of an existing transonic wind tunnel having a fixed horsepower 
(case (3)) would be about tripled if air were replaced by an argon gas 
mixture having 7 = 1.28, but only doubled by a mixture having 7 = 1 . 4 . 
The advantages of reducing 7 to the lowest practical value are apparent 
from the table. For xenon mixtures, however, there is less to gain by 
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reducing 7 below l.U. It would be desirable to conduct some experiments 
varying the proportions of a gas mixture in order to determine to what 
limit 7 can be reduced in a transonic flow and still yield data directly 
applicable to air for practical purposes . ' 


Deviations from Thermal and Caloric Perfection 


Under wind-tunnel conditions air generally is considered as being 
both thermally perfect (pm = pRT) and calorically perfect (7 = constant). 
Most other gases can not be so considered. Calculations of the degree 
of deviation from thermal perfection for the various argon gas mixtures 
are presented in Appendix A. It will suffice here to state that such 
deviations are indicated to be unimportant for normal operating conditions 
of wind tunnels. The smallness of these deviations is due in part to 
the monatomic (argon) component, which by itself is almost thermally 
perfect, and, in part, to the fact that the polyatomic component is a 
fluorochemical. The unusually low intermolecular forces in fluorochemical 
gases result in small deviations from thermal perfection. 

Deviations from caloric perfection, on the other hand, are not small 
and represent an important technical consideration. The specific heat 
of polyatomic gases varies widely, as figure 1 clearly indicates. Thus, 
a gas mixture proportioned to yield 7 = 1.4 for one Mach number will 
not yield the same value of 7 at other Mach numbers. This may limit 
the Mach number range over which a single gas mixture could be used- satis- 
factorily. Also, at a fixed Mach number, 7 will vary with temperature. 
Under wind-tunnel conditions, 7 for air is quite constant, but varies 
under flight conditions where higher temperatures are encountered. 

To aid in evaluating the importance of caloric imperfections, the 
variation of 7 with temperature has been computed for various gas 
mixtures and for air under both wind-tunnel and flight conditions. In 
these computations the spectroscopic frequencies listed in table III 
were employed to determine 7 for the gas mixtures, whereas the tables 
of reference 62 were employed for air. Instead of showing the computed 
variation as a function of some temperature parameter, an enthalpy 
parameter (H-Hjo) /(Ht-Heo) is employed which always is zero for free-stream 
conditions and always is unity for reservoir conditions, regardless of 
the nature of the gas. Some of the results are shown in figure 3. Air 
under flight conditions at a given Moo is represented by a cross - 
sectioned band because ambient temperatures in the atmosphere vary 
considerably with altitude. Gases under wind-tunnel conditions are 
represented by a single curve since a fixed reservoir temperature (about 
40° C«100° F) has been assumed. A single mixture of 25-percent CBrFg 
with 75-percent monatomic gas is considered in figure 3(a). These pro- 
portions are about right for M^o = 2.0. At M* = 2.5, 7„ is 1.44 and 
7t is 1.34. This variation does not seem excessive. At M„ = 1.3, 7 
is uniformly less than 1.4, varying between I.37 and 1.34. In view of 
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previous comments about the apparent insensitivity of transonic flows 
to much larger variations in 7, as well as the complete insensitivity 
of subsonic flows to variations in 7^ it would appear that a single mix- 
ture could be used satisfactorily from low subsonic speeds to at least 
]^ = 2 . 5 * The situation, however, is more difficult for the Mach number 
range between 2^5 Q-nd 3 * 5 ^ as figure 3 ("b) illustrates. In this case the 
single mixture is one of ^K)-percent CF4 with 60-percent monatomic gas. 

It is to be noted that the variation of 7 is not important at conditions 
close to free-stream conditions. For example, 7 does not appear within 
the body of linearized subsonic or supersonic theory, but first appears 
in second-order terms. Hence, Busemann*s second-order theory for two- 
dimensional flow has been used to estimate the range within which 7 
must be maintained in order to introduce less than 1-percent error in 
Ap/q. The boundaries of this range, indicated by dotted lines in figure 
3(b), shrink together as Moo"^ 1 and as For the moderate Mach 

numbers under consideration, it is seen that only the region of abscissa 
between about 0.2 and 1.0 is important when assessing the variations in 
specific-heat ratio. Consequently, it is deduced that the gas mixture 
approximates flight conditions at M^o = 3»5 about as well as does air 
in a wind tunnel. Such is not the case, though, for Mg^ = 2.5 (top por- 
tion of fig. 3("b))^ and it is not known whether variations near the stag- 
nation region of the magnitude shown are important. If they are, then 
the flow in an air wind tunnel also will show important differences from 
the flow in flight through air at Mach mmibers of about h or greater. This 
may be seen from figure k by comparing the values of 7 in an air wind 
tunnel (long dashed curves) with the corresponding values of 7 in flight 
(cross-sectioned region). The observed differences are considerably 
greater than any of figure 3 ('b) between flight conditions and gas-mixture 
wind-tunnel conditions. It is also evident from figure h that the tem- 
perature dependence of 7 for polyatomic gases may actually represent 
a real virtue, particularly in hypersonic research. For polyatomic gases 
under wind-tunnel conditions, as well as for air under flight conditions, 
these caloric imperfections arise from the same physical phenomenon, 
namely, temperature dependent energy of vibration between atoms within 
a molecule. Under flight conditions at Mach numbers of about 6 or greater, 
additional energy can be transferred through dissociation, which is 
responsible for values of 7 less than 1.28 in flight. 

In order to achieve hypersonic Mach numbers in a wind tunnel without 
supersaturating the test medium, it is necessary, of course, to heat the 
mediiam. The approximate reservoir temperatures (estimated to nearest 
25° K) corresponding to the particular curves in figure k, are as follows: 


Reservoir temperature ^ to prevent 
saturation; p^ = 100 psia 


Gas 

M« = 5 

Moo = 6 

M<„ = ■ 

Air 

350 

450 

550 

0 .U-CF 4 + 0 . 6 -AT 

475 

600 

675 

0 . 2 -CF 4 + 0 . 2 -AT+ 0.6 air 

500 

625 

750 
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As might be expected, the gas mixtures require higher operating temper- 
atures than does air. It is interesting, however, that this trend can 
be reversed if it is not required that the caloric imperfections of gas 
mixtures sensibly duplicate those of air under flight conditions. For 
exajnple, to achieve a Mach number of 7 at 100 psia reservoir pressure 
would require air to be heated to about 550° K, as tabulated above, but 
would require pure CF4 (7 between 1.12 and 1.33) to be heated only to 
about 4500 K, and pure CF3CF3 (7 between I.08 and I.I9) to be heated 
only to about 375 ° K. This trend is opposite to what would be expected 
from consideration solely of the respective boiling points (-l80° C for 
air, -128° C for CF4, and C for CF3CF3) because of the dominating 

effect of the reduced values of 7. Reducing 7 will reduce the rate 
at which temperature varies with pressure in an expansion process, since 

( ^-1) /y 

T»vp . This results in producing lower pressures, and hence higher 

Mach numbers, before the temperature is reduced to a point where saturation 
occurs. From a physical viewpoint, this phenomenon can be traced to the 
conversion of internal vibrational energy of a molecule into directed 
kinetic energy during expansion of a polyatomic gas. Thus, for a given 
reduction in random kinetic energy (given temperature drop), the greater 
the initial internal vibration energy the greater will be the directed 
kinetic energy at completion of expansion. Hence, to achieve a given 
ratio of directed to random energy (given Mach number), the polyatomic 
gas need not be expanded to as low a temperature. 


Deviations from Thermal Equilibrium 
(Heat Capacity Lag) 


When a polyatomic molecule passes through a region of rapid change 
in temperature, the energy stored in vibration between atoms within the 
molecules does not always adjust to its environment with sufficient rapid- 
ity to maintain thermodynamic equilibrium. This leads to a time lag in 
the heat capacity determined by the relaxation time p. Measurements 
show that 3 increases if either pressure or temperature decreases. An 
excessive long relaxation time could significantly affect a high-speed 
flow, especially in a wind tunnel where the pressures' and. temperatures 
encountered are low, and can magnify such effects. 

Detailed computations are presented in Appendix B which provide an 
estimate of relaxation-time effects for various gases. Here only end 
results are discussed. The relaxation time (Pa) at atmosphere conditions 
is known for some of the fluorochemical gases, and generally is the order 
of 10 second. For others, 3^ is not known from direct measurements, 
but can be estimated from related measurements by using certain empirical 
relationships between relaxation time and molecular structure. Thus, Pa 
for CBrFa is estimated to be in the range 10"® to 10""^ second, and that 
for CF3CF3 to be in the range 10"® to 10 second. The relaxation time 
for dry air, on the other hand, is relatively long, being the order of lO-s 
second. 
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Relaxation time is not the only important quantity which must be 
considered. The basic parameters which measure effects of heat-capacity 
lag on a gas flow involve the product PCi^ where Ci is the heat capac- 
ity of internal vibrations within the molecules. From the estimates in 
Appendix B, heat-capacity lag appears of more importance in altering 
boundary-layer flow than in altering pressure distribution, and is of 
greatest importance at low Reynolds numbers. For air under wind-tunnel 
conditions, P is relatively long, but is negligible, thereby render- 

ing relaxation effects negligible even at low Reynolds nijmbers. For most 
of the gas mixtures under wind-tunnel conditions, the estimated relaxation 
effects are negligible at Reynolds numbers of the order of 10® or greater, 
but are significant at Reynolds numbers of the order of lO'^ or less. It 
should be remembered that for air under flight conditions, Cj, is not 
negligible as it is in a wind tunnel. Thus, relaxation effects can be 
significant in fli^t, although the estimates of Appendix B suggest this 
may be the case only at high altitudes (low Reynolds numbers) and at high 
flight speeds. Heat-capacity lag is estimated to be significant at 
Reynolds numbers below about 10® for moderate supersonic speeds, and below 
about 10® for moderate hypersonic speeds. Although relaxation-time effects 
in flight would not be simulated in an air wind tunnel, it is possible 
to simulate them approximately in a gas -mixture wind tunnel by employing 
two polyatomic gases in the mixture - one having pc^ greater than for 
flight through air (e.g., CF4) , and the other having PC^ less than for 
flight through air (e.g., CF3CF3) . This possibility may be of interest 
in low-density, high Mach number research. 

It should be mentioned that the idea of employing more than one 
polyatomic gas in a mixture can be utilized profitably in other ways. 

For example, various pairs of some polyatomic gases when mixed in special 
proportions form an azeotropic mixture; that is, a mixture for which the 
boiling temperature is not intermediate between the boiling temperatures 
of the two constituent gases, as ordinarily is the case, but is lower 
than the boiling temperature of either constituent. A mixture of 32- 
percent C3F8 with 68-percent CHgFa, for example, boils at -58° C, whereas 
the respective individual gases boil at -38° C and - 52 ° C (ref. 63) . Thus, 
although neither C3FQ nor CHgFg could individually be used at =2.5 
under normal wind-tunnel operating conditions without danger of condensa- 
tion (and, hence, have not been included in table IV(b)), their azeotropic 
mixture would be usable under such conditions . Other gas pairs known 
to form azeotropic mixtures are CHF3 with C2F0, CF^Cl^ with CF2CFCF3, 

CHCIF2 with CF2CFCF3, and CF2CI2 with CH3CHF2 (see refs. k8 and 63). A 
related point worth mentioning is that certain polyatomic gases which 
would not be usable singly because of flammability characteristics may 
be usable when mixed with other gases. An extreme example illustrating 
this point is the obviously flammable butane (C4HX0) which, when mixed 
with CF2CI2 in portions up to about 30 percent, is no longer flammable 
(ref. 64) . 
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CONCLUDING REMARKS 


By employing in place of air a test medium consisting of a heavy 
monatomic gas mixed with a heavy polyatomic gas in proportions that yield 
the desired value of 1.4 for the ratio of specific heats, it appears pos- 
sible in conducting aerodynamic research to realize certain significant 
advantages. The existence of some of these advantages has been demon- 
strated by previous investigations concerned with pure, heavy, polyatomic 
gases, especially the Freons. Because pure polyatomic gases have a 
specific -heat ratio much lower than that of air, they cannot be employed 
to yield aerodynamic data directly applicable to. air for flows involving 
extensive supersonic regions. - Aerodynamic data obtained with the gas 
mixtures considered herein would directly apply to air. 

A salient result of this study is the observation that a gas mixture 
can be concocted which behaves - under low-temperature wind-tunnel condi- 
tions - dynamically similar in several respects to air under high- 
temperature fli^t conditions, whereas air under low- temperature wind- 
tunnel conditions behaves dissimilarly. This situation arises because 
the flow of two gases can be made dynamically similar on a macroscopic 
scale, even though they differ in microscopic structure, provided all the 
pertinent, dimensionless , macroscopic parameters (such as involve, for 
example, relaxation phenomena and temperature variation of specific heat) 
are duplicated between the two gases. In achieving dynamic similarity, 
it is not necessary to duplicate any dimensional quantity such as temper- 
ature . 


Many of the fluorochemical gases studied in the present research are 
more inert than Freon 12. This extreme inertness is highly desirable, 
especially since some of these gases when mixed with the necessary small 
amount of oxygen can be breathed even in large quantities without observ- 
able effect by animals, and presumably also by man. Likewise, because 
of their inertness they could be used without undesirable effect on 
machinery. 

Since the magnitude of the advantages that can be achieved through 
use of gas mixtures has been shown to be large, they may outweigh the 
evident operational disadvantages which result from use of any gas other 
than air. Whether or not the advantages surmount the disadvantages can 
only be ascertained by detailed study of individual cases. Clearly, 
though, the possible use of gas mixtures warrants consideration in the 
design of future research facilities. The' advantages of gas mixtures also 
may be worth consideration in regard to certain existing facilities where 
the Reynolds number and/or Mach number range is now undesirably low, or 
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where the annual power consumption is undesirably high. In view of the 
magnitude of the advantages, it appears that some experimental work with 
gas mixtures is in order. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. 12, 1954 
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APPENDIX A 

ESTIMATE OF DEVIATIONS FROM THERMAL PERFECTION 


By definition, a thermally perfect gas obeys the equation of state 
pm = pRT given previously as equation (2) . Real gases closely follow 
this equation if the density is low, or if the temperature is near the 
Boyle temperature (about 2.6 T^,), but deviate from it under other condi- 
tions. A more exact state equation is 


pm 

PRT 



(Al) 


where B(T) is the second virial coefficient having dimensions of volume. 
According to the law of corresponding states, B(T) is proportional to the 
critical volume V^. 

From the results of reference 65 ^ the second virial coefficient for 
a pure gas is approximately 0.75VcB*(T*) , where T* =1.3 T/Tc and where 
B*(T*) is a fixed dimensionless function tabulated in reference 65* For 
a mixture of two gases the appropriate relations are: 


B(T) = Xi%(T) + 2 xiX2Bi2(T) + X2^B2(T) (A2) 


where 


Bi = 0.75 VciB*(1.3T/Tci) 


Ba = 0.75 Vc2B*(1.3T/Tca) 


> (A3) 


Bi 2 = 0.75 




B*(1.3T/VTciTc2) 


By use of these equations, the term pB(T)/RT of equation (Al) - which 
represents the fractional deviation from thermal perfection - has been 
calculated for the gas mixtures listed in table IV. Computations were 
made both for reservoir conditions (pt>Tt) and test-section conditions 
(p^,T^) at Mach numbers of 1.3^ 2.5^ and 3*5* For Moo = 1.3 the values 
of ptB(T-t-) /RT-t and PooB(loo) /RI« were about the same. The average value 
of |pB(T)/RT| per atmosphere of reservoir pressure varied from 
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about 0 . 00 l 4 (C4F0^ C3F0, CF3CF3, CF4, and NF3 mixtures) to about 0.004 
(CHgFg mixture). Inasmuch as transonic research apparatuses commonly 
are limited to reservoir pressures of several atmospheres (often because 
of critical model loads), the resulting deviations from thermal perfection 
would be within about 1 percent and can be neglected. For Meo = 2.5 and 
= 3 . 5 > the values of 1 p»B(T;o) I are substantially lower than 
|ptB(Tt)/RTt| . Both are tabulated as follows: 




= 2.5 

M 

«. = 3.5 

Gas 

1 

. R T„ 

|PtB(Tt) 1 
RTt 

1P»B(I^) 

RT„ 

1 iPtB(Tt) 1 

_ RTt 

CBrF3 

0.0017 

0.0033 



SFe 

.0013 

.0027 



CF3CF3 

.0010 

.0018 



CCIF3 

.0016 

.0030 



CF4 

.0012 

.0019 

0.0007 

0.0023 

KF3 

.0012 

.0019 

.0007 

.0022 

CHF3 

.0017 

.0034 



Air 

.0005 

.0002 

.0002 

.0002 

De deduced 

that gas 

mixtures at test 

-section 

conditions (p„ 


rTj 

would deviate 1 percent from thermal perfection at reservoir pressures 
between about 6 and 10 atmospheres for . = 2.5^ and at about ik atmos- 


pheres for = 3 * 5 « Since practical operation pressures of wind tunnels 
are well within these limits, the flow of these gas mixtures over a typical 
model can be regarded as closely approximating that of a thermally per- 
fect gas. Deviations from thermal perfection at reservoir conditions 
(pt^Tt) are two to three times as large as at test-section conditions. 
Hence, when computing test-section Mach number or dynamic pressure from 
measured values of p-^^ and T^, the deviations near reservoir conditions 
may have to be considered for reservoir pressures of the order of 5 to 
10 atmospheres or greater. 
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APPENDIX B 

RELAXATION-TIME DATA AND ESTIMAIK OF 
EFFECT FOR VARIOUS GAS MIXTURES 


Results of measurements of relaxation time at atmospheric pressure ' 
(3a) for some of the polyatomic gases considered in this report are as 
follows: 


Gas 

Temperature, 

°K 

Pa^ 

microseconds 

Reference 

CCI2F2 

293 

0.09 

66,4 

CF4 

373 

.66 

67 

CHF3 

373 

.42 

67 

CHCIF2 

293 

.10 

66 

CH2F2 

373 

.055 

67 

CO2 

293 

7.0 

66 

NgO 

293 

1.1 

66 

Air (dry) 

300 

^ 14 . 

68 


Although measurements of 3a could be found only for those gases listed, 
certain known characteristics of relaxation time enable a rough estimate 
to be made for some pf the remaining gases; for example, molecules with 
freedom of internal rotation invariably have very. short relaxation times, 
generally less than 10"® second (see ref. for example). Also the 
smaller the lowest fundamental vibration frequency, the shorter the relax 
ation time (see refs. 66 and Sj) . These general characteristics indicate 
that fluorocarbon molecules with internal rotation (e.g., CF3CBrF2, 
CFgCClFg, CFqCFq, CHqCCIF^, CHgCHF2, ^3^0^ and but not C4F0) 

would have relaxation times (3a)’ in the range 10”® to 10”® second, and 
that the methane-like gases CBrFs and CBrFClFa would have a relaxation . 
time in the range of 10 to 10"® second. It is to be noted that relax- 
ation times listed are for pure gases, and that the presence of certain 
impurities can greatly reduce 3a for those gases having relatively long 
relaxation times. 

Inasmuch as the relaxation time depends on temperature and pressure, 
some method of estimating 3 for the particular conditions encountered 
in wind tunnels must be employed. The analysis of Bethe and Teller 

^This value corresponds to measurements of 3a for oxygen in the 
presence of nitrogen but in the absence of water vapor. Dry air is 
considered since air in supersonic wind tunnels and in the stratosphere 
is essentially dry. 
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(ref. 68 ) yields for the relaxation time of the lowest frequency 
(v) mode 


0 _ ^10 

Cp - Ci '' 2(1 - 


(Bl) 


where Ci is the total contribution to specific heat of internal vibra- 
tion energy, Z is the total number of collisions one molecule experiences 
per second, and is the average number of collisions required to 

deactivate the lowest mode from the first quantum state to the zero state. 
Subsequent calculations are based on the additional equations 


p 

Z = 1.3 jr 

(ref. 70 ) 

(B2) 

3~ n 

(ref. 71 ) 

(B3) 


^ 10 - 25 (^ 1 ^)''''^ (ref. 68 ) (B4) 


where b is a constant for each gas evaluated from measurements of' 3 a 
at the temperatures listed in the above table. These equations were 
developed for pure polyatomic gases. They will be used also for gas 
mixtures since the basic theory (especially eq. (b 4)) is highly approxi- 
mate, and since the effectiveness of an argon molecule in exciting vibra- 
tion generally is roughly the same as that of an inert polyatomic molecule 
(See the summary table in ref. 72 and the specific data for CHClFa-Ar 
mixtures in ref. 73 •) 

Relaxation times 3 oo corresponding to free-stream conditions (Too,p») 
in a wind tunnel have been calculated for assxamed reservoir conditions 
of one atmosphere total pressure and C total temperature. The results 
together with the "relaxation distance" d='3»Uoo, were computed to two 
significant figures and are tabulated to one significant figure as follows 
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Gas 

« 

M„=1.3, T«=230° K, 
P„=0.36 atm 

Moo 2 . 5 , Too=l^° K, 
Poo = 0.059 atm 

Moo 3.5> K, 

Poo = 0.013 atm 

sec 

d, 

in. 

3oo> 

sec 

d; 

in. 

3oo? 

sec 

d, 

in. 

CCl^Fg 

2x10-7 

0.002 





CF4 

3X10"® 

.03 

5x10"® 

0.8 

3x10"^ 

6 

CHF3 

2X10"® 

.02 

3x10-® 

.5 



CHCIF2 

2x10-7 

.003 





CH2F2 

1X10“ 

.001 

. 7x10-7 

1 — 1 

0 

• 



CO2 

1^X10"^ 

•5 





N2O 

6x10"® 

.06 





Air (dry) 

7x10-® 

1. 

3x10"® 

60 . 

4x10"® 

1000 


It is interesting that under wind-tunnel conditions, air, has the longest 
relaxation time. For M(o>2.5 any model less than 5 feet long in an air 
wind tunnel would he entirely immersed within the relaxation region from 
the bow wave. This does not mean, however, that relaxation effects are 
important for air in a wind tunnel because the amount of Internal vibra- 
tion heat capacity (Ci) that lags is very small. The effect of heat- 
capacity lag on flow conditions is influenced by both Cj^ and 3* 

To determine the dimensionless parameters which form a pertinent 
measure of relaxation effects, use is made of the theoretical equivalence 
between relaxation phenomena and volume-viscosity phenomena. In refer- 
ences 'JO, 'jk, and 75> this equivalence is shown to be valid at least under 
conditions where the relaxation time is small compared to the character- 
istic time (l/uoo for aerodynamic problems), and where Ci is small com- 
pared to Cp. For purposes of estimation, then, the complete, first-order 
stress tensor, including the volume viscosity, is used as the basis of 
computation . 

8uj^ 

PlJ = (I'lj'Mavier + ' (®5)' 

Stokes 


In this equation Hy coefficient of volume viscosity, 8uj,./8xjj. 

is (using summation convention) the divergence of the velocity vector, 
and 6ij is unity if i is equal to j, but is zero otherwise. The 
ultrasonic theory of Tisza (ref. 7^)> the kinetic theory of Kohler 
(ref. 70), and the irreversible thermodynamics theory of Meixner (ref. 75) 
each indicate Py to be proportional to 7(7-l)p3^ although the propor- 
tionality factor f = My/7(7-l)p3 differs somewhat between these analyses. 
The various values calculated for f are Ci/Cp, Ci/(Cp- 7 Ci), and 
Ci/(Cp-7Ci), respectively. These values are mutually consistent in view 
of the common assumption made in references 'Jk, "JO, and 75> that is 
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dmall compared to Cp. For present purposes an intermediate value 
f= Cp/(Cp-Ci) is used^ since a comparison with the Kantrowitz impact- 
pressure theory (ref. 76) indicated this value to yield the same energy 
dissipation - and, hence, the same total -pressure defect - as the 
Kantrowitz theory for a gas flow wherein small relaxation effects occur. 
Consequently, the evaluation Py = 7(7-l)pPCi/(Cp-Ci) together with equa- 
tion (B5) will he employed to estimate effects of heat-capacity lag both 
for inviscid and viscous flow. 

For inviscid flow the Wavier-Stokes portion of the stress tensor is 
simply the static pressure p. The pertinent similarity parameter involv- 
ing relaxation time is hyuj'(pl) = 7(7-1) (3ua,A)Ci/(Cp-Ci) . Since the 
equation of steady motion is pujBuiAxj = the ratio 

of the relaxation term to the other stress term is, in order of magnitude. 



(B6) 


Which is seen to be the same as the similarity parameter.^ Hence, 
is interpreted as a measure of the fractional effect of relaxation phenom- 
ena on pressure distribution. 

^ Inasmuch as the parameter iiyU<„/pZ can be rearranged as (pv/ii) 

7M«o /He, it is clear that relaxation effects on pressure distribution 
generally will not be important except for low-density, high Mach number 
conditions where Re is small. For dry air at wind-tunnel temperatures, 
My/li is negligible compared to unity, and hence the parameter |i u„/pZ 
also is negligible, even at low Reynolds numbers. For dry air at^f light 
temperatures (say 300° K to 2000° K), however, varies between about 

80 and 200, thereby indicating sizable relaxation effects, for example, 
at = 3*5 if the Reynolds number of flight is the order of several 
thousand or less. For gas mixtures at wind-tunnel temperatures, p^/p 
varies widely depending on the polyatomic gas , and covers the range from 
about 20 to 1000, thus encompassing values both less than and greater 
than those for air under flight conditions . 

In a gas mixture, Ci = xaCia. Computations of the resulting para- 
meter p^u'^p^gZ = 7(7~l) (Poo^/Z)Ci/(Cp-Ci) for wind-tunnel flow conditions 

A more precise calculation of this ratio of stress terms can be made 
using the continuity equation to eliminate ^ukAxk. This procedure ulti- 
mately yields a value (HyU„/pZ)[M2a(u/u«,)/a(x/z) ], from which it can be 
deduced that, the simple parameter (PyU^/pZ) will underestimate relaxation 
effects for bluff obstacles like the sphere, but will overestimate them 
for slender obstacles like an airfoil. 
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over a 1-inch model yield the following results: 


Gas 

r— — 

/ V XoC ^ 00^^00 ' 

Cp-C- 1 ^ 1 Pt - 1 

Woo = 1-3 

m„ = 2.5 

LTN 

. 

on 

II 

CC 1^2 

0.0003 



CF4 

0 

0 

• 

o.oU 

0.1 

CHF3 

.002 

.01 


CHCIF2 

.0003 




.00006 

.0001 


CO2 

.03 



N2O 

.005 



Air (dry) 

.0002 

. ■ M.i - n— ,1. 

.00005 

.00008 


For transonic wind tunnels it is seen that only CO 2 mixtures would exhibit 
appreciable relaxation effects on pressure distribution. On the other 
hand, for mixtures involving CF^ would be affected roughly by 

10 percent for the assumed conditions of p-j^ = 1 atmosphere and 2=1 inch. 
To maintain these effects below about 1- percent would require the product 
p^2 to be greater than 10 atmosphere-inches (e.g., a 5~inch model in a 
2-atmosphere tunnel) • Such restrictions are believed not to represent 
serious limitations on the usefulness of gas mixtures, especially since 
the favorable effect of certain impurities in reducing p has not been 
included in the analysis. It is concluded, then, that insofar as pres- 
sure distribution in inviscid flow is concerned, relaxation effects need 
not be considered for the usual wind-tunnel operation conditions. 

For viscous flow within a laminar boundary layer, the Navier-Stokes 

portion of the stress tensor customarily is approximated by . The 

dy 

ratio of the relaxation term to the ordinary term in the equation of 
motion is 



Examination of the solution to boundary-layer equations for flow over a 
flat plate (e.g., ref. 77) indicates the average value of (Suj^/Sxj^) /(Su/^y) 
to be of the order of (7-1 )Moq^/iO ^/Re, and b/l to be of the order of 
[6 + Hence, 
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This ratio, which is interpreted as a measure of the fractional effect 
of relaxation phenomena on houndary-layer flow, is important primarily 
at low Reynolds numbers and high Mach numbers, just as is the effect of 
relaxation phenomena on inviscid flow. For air under wind-tunnel condi- 
tions, the expression (B7) is entirely negligible due to the very low 
values of However,, for air under flight conditions (dv/|i~10^), 

the ratio can be significant at low Reynolds numbers (e.g., at Re~500 
for Meo = 3.5^ at Re~5000 for M« = 7^ and at Re ''-200,000 for 
Moo = 10) . In order to simulate in a low-density gas-mixture wind tunnel 
the heat-capacity lag that would exist in the boundary layer during flight, 
the ratio would have to be simulated. As previously noted, it is 

possible to simulate this ratio in a gas-mixture wind tunnel due to the 
wide range of values of iiy/n covered by the various polyatomic gases. 
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a 

c 

Cp 

CR 

g 

h 

R 

HP 

k 

I 

L 

m 

M 

n 

P 

<1 

R 

Re 

T 

u 

V 

w 


APPEHDIX C 
NOTATION 

speed of sound 

speed of light 

specific heat per unit mass 

specific heat per mole 

wind-tunnel compression ratio 

acceleration of gravity 

Plank’s constant 

molar enthalpy 

horsepower 

Boltzman's constant 

characteristic length of model 

characteristic dimension of wind-tunnel test section 
molecular weight (29 grams per mole for air) 

Mach number 

number of atoms per molecule 
pressure 

1 5- 

dynamic pressure, — Pu*^ 

universal gas constant per mole (I.987 cal. mole"^ °K"^) 

Reynolds number 

temperatiore 

velocity 

volume of gas 

mass-flow rate, (pu) times (cross-section area) 



. mole fraction, equal to ratio of partial pressure to mixture 
static pressure and equal to fraction by volume 

relaxation time 

ratio of specific heats 

mass density 

viscosity coefficient 

volume viscosity coefficient 

wave number (cm”^, from spectroscopic data) 

Subscripts 

free-stream conditions in test section 

critical conditions for gaseous phase of a compound 

total conditions for gas brought isentropically to rest (wind- 
tunnel reservoir conditions) 

0° centigreide 

monatomic gas 

polyatomic gas 

Superscripts 

quantity divided by corresponding quantity for air 
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TABLE I.- PROPERTIES OF HEAVY MONATOMIC GASES AND AIR CONSTITUENTS 


Symbol 

Gas 

m, 

gm.. 

mole 

Boiling 

temp., 

°C 

Critical 

constants 

Viscosity 
at 0 ° C 
^^o^ 

micropoise 

To 

°K 

Pc^ 

atm. 

cc/gm. mole 

Ar 

argon 

k) 

-186 

151 

48 

75.3 

210 

Kr 

krypton 

84 

-153 

210 

54 

92.1 

230 

Xe 

xenon 

131 

-107 

290 

58.2 

ll 4 

210 

Ns 

nitrogen 

28 

-196 

126 

33.5 

90.0 

166 

O2 

oxygen 

32 

-183 

154 

49.7 

74.4 

189 


air 

29 


132 

37.2 

83 

172 




TABLE II.- GENERAL CHARACTERISTICS OF POLYATOMIC GASES 


R 
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hi 





TABLE III.- PHYSICAL CHARACTERISTICS OF POLYATOMIC GASES 
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k2 




estimated from similar gases, or from liquid density and reduced orthobaric density of similar gases- 
number in parentheses denotes degree of degeneracy of wave number it follows. 




TABLE IV.- CHARACTERISTICS OF GAS-MIXTURE WIND OUNNELS 
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Figure 1.- Some typical examples of specific heat computed from spectroscopic 
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Figxire 2.- Some typical examples of viscosity coefficient for gas 

mixtures. 
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Figure 4 .- Variation of specific-heat ratio between free-stream and 
stagnation, conditions at hypersonic speeds. 
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